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I.  INTRODUCTION 


"  The  purpose  of  this  document  is  to  define  and  briefly  discuss  coordinate 
systems  commonly  employed  at  test  ranges  as  well  as  to  recommend  procedures  and 
practices  for  the  use  of  a  global  coordinate  systemic  This  document  replaces 
IRIG  Document  151-69  (formerly  103-69),  Global  Coordinate  .System. 

This  document  is  primarily  intended  for  use  by, individual s  familiar  with 
the  basic  concepts  and  definitions  commonly  employed  in  geodesy.  In  order1 to 
aid  the  reader  who  requires  additional  background  information,  a  bibliography 
of  some  of  the  standard  texts  on  geodesy,  map  projections,  and  spherical 
astronomy  is  provided.  — ^  (  ^  *  ....  X  ■  -'*  i 

This  document  specializes  in  those  trajectory  coordinate  systems  that  are, 
or  could  be,  potentially  applicable  to  data  reduction  tasks  performed  at  the 
various  test  ranges.  No  attempt  has  been  made  to  address  the  subject  of  coor¬ 
dinate  systems  in  a  general  manner  or  to  develop  a  compendium  of  the  various 
coordinate  systems  employed  at  the  test  ranges.  General  discussions  on  coordi¬ 
nate  systems  are  given  in  a  number  of  standard  texts  on  geodesy  and  are  not 
repeated  in  this  document.  Information  on  special-purpose  coordinate  systems 
employed  at  a  given  range  can  usually  be  obtained  from  documents  at  the  local 
installation.  A  compendium  of  all  coordinate  systems  used  by  member  ranges 
would  produce  a  voluminous  document  of  dubious  value.  In  spite  of  being 
narrow  in  scope,  this  document  still  manages  to  serve  a  useful  purpose 
through  an  in-depth  presentation  of  the  definitions,  properties  and  useful 
transformations  associated  with  those  coordinate  systems  having  a  general 
interest  to  data  reduction  specialists  at  most  test  ranges. 


II.  GLOBAL  COORDINATES 

Any  global  coordinate  system  (geodetic  datum)  is  based  on  an  ellipsoid  of 
revolution  whose  size  and  shape  approximate  that  of  the  entire  Earth  and  whose 
center  coincides  in  some  manner  with  the  gravitational  center  of  the  Earth. 
Examples  of  such  systems  are  the  Mercury  Datum  based  on  the  Fischer  1960 
Ellipsoid  and  the  DOD  World  Geodetic  System  (WSS)  1972. 

Satellite  and  missile  trajectories  are  relativa  to  the  Earth's  center  of 
gravity;  therefore,  range  data  referred  to  toe  center  of  gravity  are  the  end 
products  for  such  trajectory  analyses.  When  the  sensor  coverage  exceeds  the 
extent  of  local  or  regional  datums,  the  sires  must  be  located  on  a  global  datum 
to  achieve  consistent  geometry. 

Global  coordinates  based  on  a  given  ellipsoid  are  determined  by  transforma¬ 
tions  from  regional  or  local  coordinate  systems  or  from  observations  of  a  satel¬ 
lite  whose  orbit  is  known  relative  to  seme  coordinate  system:  These  computations 
and  transformations  are  usually  assumed  to  lead  to  consistent  results.  However, 
when  it  is  possible  to  make  comparisons,  the  results  are  rarely  in  complete 
agreement.  Differences  in  computational  and  orbital  force  field  models,  along 
with  parameter  and  computer  characteristics,  are  responsible  for  such  discrep¬ 
ancies.  Additionally,  in  theory,  the  centers  of  two  global  datum  systems 
should  only  differ  by  a  constant;  hewever,  in  practice,  the  difference 
appears  to  be  a  function  of  the  location  of  the  measurements. 


III.  REGIONAL  COORDINATES 


A  regional  geodetic  datum  consists  of  an  ellipsoid  whose  size  and  shape 
approximate  that  of  the  particular  region  of  the  Earth.  This  ellipsoid  is 
oriented  to  a  physical  point  in  the  region,  and  this  point  is  the  origin  of 
the  regional  datum.  Examples  of  such  systems  and  their  reference  ellipsoids 
are  the  European  Datum  (ED)  on  the  International  Ellipsoid  of  1924,  the  North 
American  Datum  (NAD)  on  the  Clarke  1866  Ellipsoid,  and  the  Tokyo  Datum  (TD)  on 
Bessel  1841  Ellipsoid. 

Conventional  surveys  in  a  region  are  referenced  to  the  datum  network  defin¬ 
ing  the  regional  system.  When  distance  and  azimuth  between  regional  sensors 
are  necessary  for  seal ing,  such  regional  coordinates  are  the  best  computational 
source.  Regional  coordinates  may  be  extended  to  greater  distances  by  optical 
and  electronic  observation  of  a  satellite.  The  relation  between  regional  and 
global  coordinate  systems  is  specified  by  displacement  between  ellipsoid  centers 
in  rectangular  coordinates  (AE,  AF,  AG)  and  differences  in  ellipsoid  sizes  and 
shapes  (Aa,  Af),  where  a  and  f  are  semi-major  axis  and  flattening,  respectively. 


IV.  PARTICULAR  RANGE  (LOCAL)  COORDINATES 

These  coordinate  systems  are  usually  rectangular  systems  whose  center 
(origin)  and  x-y  plane  are  tangent  to  some  convenient  point  on  the  range 
(usually  one  of  the  launch  pads  or  instrumentation  sites).  The  coordinates 
may  be  transformed  from  these  local  coordinates  to  appropriate  regional  coor¬ 
dinates  or  global  systems.  These  local  range  systems  are  primarily  a  compu¬ 
tational  convenience  designed  for  the  range  configuration'  and  mission.  Such 
systems  are  not  amenable  to  interrange  projects  or  missions. 


V.  DATUM  TRANS  FCSMATI'ONS 


Frequently,  it  is  neces-'ry  to  transform  station  geodetic  coordinates  from 
one  datum  to  another.  The  following  information-  is  necessary  for  this  transfor¬ 
mation  : 


a0  -  the  semi -major  -xis  of  the  reference  ellipsoid  of  the  datul^^m  which 
the  transformation  will  occur  (original  datum), 

fo  =  the  flattening  of  the  reference  ellipsoid  of  the  original  datum, 

a^  =  the  semi -major  axis  of  the  reference  ellipsoid  for  the  datum  to  which 
the  transformation  will  occur  (new  datum), 

f^  =  the  flattening  of  the  reference  ellipsoid  of  the  new  datum,  and 


AE  =  En  -  E0 
AF  =  Fn  -  F o 

AG  =  Gn  -  Gq 


the  displacement  of  the  center 
of  the  new  datum  relative  to  the 
original 


The  EFG  system  is  a  right-handed  Cartesian  system.  The  origin  of  this 
system  is  at  the  geometric  center  of  the  reference  ellipsoid.  The  coordinate 
axes  areorionted  as  follows:  E  and  F  ,1 ie  in  the  plane  of  the  Equator,  G 
coincides  with  the  rotational  axis  of  the  Earth  ard  is  positive  througn  the 
North  Pole,  E  is  positive  through  Greenwich  Meridian  (0°  longitude),  and  F  is 
positive  to  complete  a  right-handed  system.  This  system  rotates  with  the  Earth 
and  is  Earth-centered  and  Earth-fixed  (ECEF).  The  fact  that  the  various  Earth- 
centered  systems  exhibit  displacement  of  their  centers  is  evidence  of  experi¬ 
mental  error  rather  than  disagreement  in  fundamental  assumptions.  This  EFG 
system  will  he  called  the  Master  Coordinate  System  or  the  system  in  which 
integration  takes  place. 

The  method  of  transformation  consists  of  the  following  steps: 


1.  Calculate  the  geocentric  Cartesian  coordinates  of  the  original 
coordinates  given  4>q, '  Xg,  H0,  a0,  f0,  where  <fr0,  Xg  and  H0  are  the 
geodetic  latitude,  longitude  and  height  above  the  reference 
ellipsoid,  respectively,  of  the  point  to  be  transformed.  The 
method  of  transformation,  discussed  in  section  VII,  paragraph  7, 
produces  the  data  set  (E0,  F0,  G0). 

2.  Calculate  the  geocentric  Cartesian  coordinates  of  the  hew  datum: 


En  *  E0  +  a£ 
»  F0  ♦  AF 
Sfl  *  So  +  AG 


\  IS  BIANK 
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3.  Calculate  <f>^,  X^,  and  H.^  using  a^,  f^,  E^,  F^,  and  as  shown  in 
section  VII,  paragraph  7. 


The  origins  and  ellipsoids  for  a  few  datums  are 


NAD  1927 

ED 

TO 

Meade's  Ranch,  KS 

Potsdam 

Tokyo 

4> 

39°13'26.6860 "N 

52°22,54.81"N 

35°39' 1 7.51 "N 

X 

261 °27'29. 49400" 

13°04‘01 .66" 

139°44,40.50" 

H 

599.4M±0.3 

Ell  ipso id 

Clarke  1866 

International  1 924 

Bessel  1841 

The  ellipsoid  parameters  of  interest 

are 

Ellipsoid 

a  (meters) 

f 

Clarke 

1866 

6378206.4 

1/294.9786932 

Fischer  (Mercury)  1960 

6378166 

1/298.3 

Kaula 

1  961 

6378165 

1/298.3 

6377397.155 

6378135 


Bessel  1841 
WGS  72 

International  (Hayford)  1924 


63^8388 


1/299.1528128 

.1/298.26 

1/297.0 


VI .  TRAJECTORY  COORDINATE  SYSTEMS 


Orbital  computations  involve  a  number  of  different  coordinate  systems 
and  transformation  of  data  among  these  systems.  These  coordinate  systems  are 
related  to  each  other  through  a  reference  figure  representing  the  actual  Earth. 
This  reference  figure  is  an  ellipsoid  of  revolution  of  given  dimensions  whose 
surface  is  assumed  to  approximate  closely  the  mean  sea  level  surface  of  the  . 
actual  Earth.  It  is  also  assumed  that  the  inass  of  the  Earth  is  sufficiently 
homogeneous  for  its  center  to  be  taken  at  the  geometric  center  of  the  reference 
ellipsoid.  In  this  way,  the  major  axis  of  the  reference  ellipsoid  lies  in  tne 
equatorial  plane  of  the  Earth,  and  the  minor  axis  coincides  with  the  Earth's 
rotational  axis. 

The  location  and  origin  of  a  given  data  acquisition  coordinate  system  are 
determined  by  an  astronomical  and  a  geodetic  survey  (the  latter  of  which  is 
based  upon  a  specific  ellipsoid).  The  geodetic  coordinates  a^e  latitude  <p, 
longitude  X,  and  height  H  above  the  surface  of  the  reference  ellipsoid.  The 
geodetic  latitude  of  a  point,  which  is  positive  in  the  Northern  Hemisphere, 
is  the  angle  between  the  equatorial  plane  and  the  geodetic  vertical  through 
the  point  extended  to  intersect  the  plane  of  the  equator.  The  geodetic  longi¬ 
tude  of  a  point  is  an  angle  in  the  equatorial  plane  with  its  vertex  at  the 
center, of  the  Earth,  its  initial  side  through  the  Greenv;ich  Meridian,  and  its 
terminal  side  through  the  meridian  of  the  point.  Longitude  is  measured  posi¬ 
tive  eastward  of  Greenwich.  Geodetic  height  is  measured  positive  up  from  the 
surface  of  the  reference  ellipsoid  along  the  geodetic,  vertical . 

The  astronomical  coordinates  are  latitude  and  longitude  X'.  The  astro¬ 
nomical  latitude  is  defined  as  the  angle  between  the  equatorial  plane  and  the 
astronomical  vertical  through  the  point.  The  astronomical  vertical  is. defined 
by  the  direction  of  gravity  (plumbline)  at  the  point.  The  astronomical  longi¬ 
tude  is  the  angle  between  the  plane  containing  the  Greenwich  Meridian  and  the 
projection  of  the  astronomical  vertical  on  the  equatorial  plsne.  The  latitude 
and  longitude  differences  between  geodetic  ar.d  astronomical  coordinates  define 
the  deflection  of  the  vertical.  Astronomical  coordinates  define  only  orienta¬ 
tion,  while  geodetic  coordinates  define  orientation  and  location.  This  dis¬ 
tinction  is  of  importance  to  the  transformation  equations  detailed  in  the  next 
section. 

Prior  to  any  orbital  computation,  a  set  of  initial  conditions  must  be 
obtained  if  equations  of  motion, are  to  be  integrated.  Position  and  velocity, 
are  usually  needed  (and  ballistic  coefficient  and  covariance  matrix,  if  pos¬ 
sible).  The  initial  conditions  may  be  transformed  to  the  ECEE  system. 


VII.  PARTICULAR  TRAJECTORY  COORDINATE  SYSTEMS 

•  •  •  . 

1.  Topocentric  Range  Coordinate  System  (x,  y,  z,  x,  y,  z) 

A  topocentric  Cartesian  system  is  a:>sumed  with  the  x-y  plane  normal 
to  the  geodetic  normal  determined  by  the  origin,  with  geodetic  coor¬ 
dinates  $,  X  and  H.  The  xyz  range  system  is  shown  in  figure  1-1. 

The  coordinate  system  is  right-handed  with  positive'  x  pointing 
east,  positive  y  pointing  north,  and  positive  z  pointing  up  along 
the  normal  completing  the  right-handed  system. 


components  are 

x  =  R  cos  e  sin  a 

y  =  R  cos  c  cos  a  ,  0  ) 

z  =  R  sin  e 


where 

R  is  the  slant  range  from  the  origin  to  a  point  with  topocentric 
coordinates  (x,y,z) , 

is  the  elevation  angle  of  the  position  vector  above  the  x-y 

plane,  and 


a  is  the  azimuth  angle  measured  in  the  x-y  topocentric  plane  from 
True  North  to  the  projection  of  the  position  vector  on  the  x-y  plane. 


Hat 

—  **  blank 
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If  the  components  of  the  position  vector  x,  y,  z  are  known,  then  R,  ci,  i 
may  be  determined  from  the  following  equations- 


R  =  (x2  *  y2  +  z2)  2 


12) 


(Double  Arc  tan) 

(If  =  a<0,  then  a  =  a  +  2-r) 


(3) 


£  3  tan 


( x  ‘ 


(4) 


In  general,  position  and  range  data  are  the  basic  measurements  received  oy 
instrumentation  systems  that  are  external  to  a  vehicle  in  flight.  Most  methods 
of  obtaining  velocity  and  acceleration  make  use  of  numerically  differentiated 
filtered  position  data  (expressed  as  a  function  of  time).  If  the  position  and 
velocity  components  are  known,  the  range  rate,  azimuth  rate,  and  elevation  rate 
may  be  computed  as  follows: 


XX  +  yy  +  zz 

(5) 

R 

yx  -  *y 
x2  >  y2 

(radians 

per  second) 

(6) 

(x2  ♦  y2)4 

(radians 

per  second) 

(7) 

2.  Geocentric  Coordinate  System  (Master)  (E,  F,  G,  E,-  F,  G) 


We  assume  a  right-handed  Cartesian  system  located  at  the  center  of  the 
ellipsoid  of  reference  with  the  E-F  plane  fixed  in  the  equatorial  plane  (that 
is  ECEF)  and  the  G  axis  along  the  rotational  axis  of  the  Earth;  E  is  positive 
in  the  direction  of  longitude  (A)«  0®  and  latitude  (o)  ■  0°.  F  is  positive 
90°  counterclockwise  from  E,  and  G  is  positive  in  the  direction  $  *  90°  (see 
figures  2-1  and  2-2). 
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Prior  to  any  transformation  of  topocentric  coordinates,  corrections  for  deflec¬ 
tion  from  the  vertical  must  be  applied.  Assuming  these  corrections  have  been 
made,  the  transformation  of  position  and  velocity  components  measured  in  the 
topocentric  range  system  is  given  by 


E 

-sin 

X 

0 

-cos 

X 

0 

0 

1 

0 

0 

X 

. 

E 

0 

F 

= 

cos 

Xo 

-sin 

X 

0 

0 

0 

sin 

-cos  *0 

y 

♦ 

r 

'  0 

G 

0 

0 

1 

0 

cos 

sin 

2 

G 

0 

-sin  X  -cos  X  sin  cos  X  cos  $> 

O'  o  o  o  o 

cos  X  -sin  X  sin  p  sin  X  cos  $ 

o  oo  o  o 


COS  A 


sin  a 


y  +  f0  {^) 


The  geocentric  coordinates  of  the  topocentric  reference  site  are  computed 


I 

J 


Eo  ’  (VHo)coS  cos  Xo 
Fo  *  {VHo)cOS  *o  Sin  Xo 


Go  »  [M0(1-e*)  ♦  Hq]  sin 

The  parameter  N  in  equation  (10)  is  the  transverse  radius  of  curvature 


given  by 


( 1  -e2sin2d>0) 


where 


aQ  is  the  semi-major  axis  of  the  reference  ellipsoid, 

e  is  the  eccentricity  of  the  reference  ellipsoid, 

is  the  geodetic  latitude  of  the  reference  site,  positive  north 
0  of  the  Equator,  and 

XQ  is  the  geodetic  longitude  of  the  reference  site,  positive  east. 


In  addition,  Hq  is  the  geodetic  height  of  the  reference  site  and 

h  =  n  +  h 
o 

where  n  is  the  geoid-ell ipsoid  separation  and  h  is  the  height  of  the  reference 
site  above  (or  below)  mean  sea  level. 

The  geocentric  velocities  are  found  from  the  following  transformations: 


where  K  is  the  matrix  in  equation  (12). 

If  the  geocentric  coordinates  of  the  point  are  known,  the  topocentric  coor 
dinates  may  be  found  by 


If  the  topocentric  coordinates  are  positioned  so  that  the  x  axis 
angle  from  north  other  than  90°,  the  <  matrix  becomes 


where 


K1 1  K1 2  K1 3  i 


K21  K22  K23 
K31  K32  K33 


<11  =  -sin  X  sin  a  -  cos  a  cos  X  sin 

<12  3  sin  X  cos  a  -  sin  a  cos  X  sin 

<13  3  cos  X  cos  <fr, 

<2i  3  cos  X  sin  a  -  cos  a  sin  X  sin  <fr, 

<22  3  -cos  *  cos  a  -  sin  a  sin  a  sin 

<23  3  sin  X  cos  $, 

<31  3  cos  $  cos  a, 

<32  3  cos  4>  sin  a,  and 
<33  *  sin  ♦. 

where  a  is  the  angle  from  north  to  the  position  of  the  positive  x  axis 
the  topocentric  system  (see  figure  2-3). 
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Figure  2-3 


3.  Earth-Fixed  Spherical  Coordinate  System  (r,  X,  6,  vg,  Ag,  Yg\ 

A  spherical  system  is  assumed  with  the  same  origin  as  the  ECEF  (EFG)  coor¬ 
dinate  system  (see  figure  3-1).  In  addition,  a  UVW  coordinate  system  is  also 
introduced.  The  latter  is  a  right-handed  coordinate  system  with  1)  pointing  ea't, 
V  pointing  north,  and  W  pointing  along  the  projection  of  the  geocentric  radius  r. 
The  U-V  plane  is  normal  to  the  geocentric  radius  (see  figure  3-1)'. 


Greenwich 


Figure  3-1 

The  following  parameters  can  now  be  defined  in  terms  of  these  coordinate 
systems: 

r  is  the  geocentric  range  from  the  center  of  the  reference  spheroid  to 
the  point  (origin  of  the  (UVW)  coordinate  system). 

6  is  the  declination,  positive  north.  The  declination  is  the  angle 
between  the  equatorial  plane  and  the  geocentric  range  vector. 

X  is  the  longitude,  positive  east. 

vg  is  the  Earth-fixed  total  velocity. 

A  is  the  azimuth  of  the  velocity  vector  projected  on  the  U-V  plane. 

e  The  angle  is  measured  from  north  to  the  projection  of  the  velocity 
vector.  Ag  is  commonly  called  the  heading. 

Yg  is  the  elevation  of  the  velocity  vector  above  the  U-V  plane,  y  is  com¬ 
monly  called  the  flight  path  angle.  " 


If  the  position  vector  components  of  a  point  in  space  are  given  in  the  ECEF 
system,  then  r,  X  and  6  may  be  found  from 

r  =  (E2  +  F2  +  g2)1*  (17) 

wnere  E,  F,  and  G  are  the  geocentric  coordinates  of  the  point. 


o  =  tan' 


(E2  +  F2)7  ’ 


X  =  ta 


-i  /F\  ,  (Double  Arctan) 

n  \  E /'  Uf  X<0,  then  +  2tt) 


If  the  coordinates  of  the  point  in  space  are  given  in  polar  coordinates  r, 
X,  6,  the  geocentric  coordinates  are 

E  =  r  cos  6  cos  X\ 


F  =  r  cos  6  sin  X 


G  =  r  sin  6 


If  the  geocentric  velocities  are  known,  they  can  be  transformed  from  the  EFG 
coordinate  system  to  the  UVW  coordinate  system  with  the  following  rotations: 

0,  =  X 


e2  =  90° 
e3  3  (9,0°  -  5) 


The  transformation  is 


v  =  r3r2r1  F 


where 


i 

0 

0 

COS 

02 

sin 

,92 

0 

cos 

61 

sin 

el 

0 

R3R2R1  = 

0 

cos  9- 

sin 

93 

-sin 

62 

COS 

CD 

0 

-sin 

9i 

cos 

91 

0 

0 

-sin  03 

cos 

03 

0 

0 

1 

0 

0 

l 

Substituting  the  values  of  9,,  0?,  and  0-  into  the  R-R?R,  matrix  and  matrix  mul 
ti pi) ing  gives 


8 


-sin  X 

R3R2R1  =  "sin  5  C0S  * 


cos  X 

-sin  5  sin  X. 
cos  8  sin  X 


0 

cos  6 
sin  6 


I  cos  5  cos  X  cos  8  sin  X  sin  6  | 

The  result  of  the  transformation  using  equations  (21)  and  (23)  is 

U  =  -E  sin  X  +  F  cos  X 

•  •  •  •  i  . 

V  =  -E  sin  6  cos  X  -  F  sin  6  sin  X  +  G  cos  8 

W  =  E  cos  5  cos  X  +  F  cos  5  sin  X  •'  o  sin  6 
The, transpose  of  the  R3R2R^  matrix  is 

I  -sin  X  -sin  6  cos  X  cos  6  cos  X  I 


( R  3  R  2  R 1  ^ 


cos  X  -sin  6  sin  X  cos  5  sin  X 


cos  6 


sin  <5 


u 

A 

\£ 

* 


f  -  (r3r2r1)t 


*  •  • 

If  the  geocentric  velocities  (E,  F,  6)  of  the  point  in  space  are  known, 

v  ,  A  ,  X  may  be  found  from  the  velocity  components  transformed  to  the  UVW 
e  e  e 

coordinate  system  by  the  use  of  equation  (24)  and  the  following  relations: 


•  O  .O  .0  H 

v  =  (IT  +  r  +  W) 


,  .  (Double  Arctan) 

A„  =  tan-1  (t)  (If  Aa  <0,  then  Ag  *  Ag  +  2tt) 


y  =  tan"1  (  — - , 

e  ‘  (U2  4  v2)\ 


(27) 

(28) 

(29) 


If,  on  the  other  hand,  the  total  velocity,  heading,  and  flight  path  angle 
are  known  in  the  UVW  coordinate  system,  then 

U  =  vg  cos  Ye  sin  Ag  (30) 

V  *  v^  cos  y„  cos  Aa 

e  e  e 

W  *  ve  sin  Yg 

The  geometric  velocities  may  be  obtained  in  this  case  by  transforming  the 
U,  V,  W  to  the  E,  F,  G  (ECEF)  system  with  equation  (26).  . 

•  •  • 

4.  Cartesian  Inertial  Coordinate  System  (Ej,  Fj ,  Gj,  Ej,  Fj,  Gj) 

In  figure  4-1,  an  Earth-centered  inertial  (ECI)  coordinate  system  is  repre¬ 
sented  by  Ej,  Fj,  Gj.  The  ECEF  rotating  coordinate  system  is  represented  by 
E,  F,  G.  The  ECEF  system  is  rotating  about  the  Gj(G)  axis  with  an  angular 
velocity  9.  (9  =  uiQ,  the  rotational  rate  of  the  Earth.)  By  definition,  Ej  is 

positive  toward  the  mean  vernal  equinox  of  date,  Fj  is  positive  90*  counter¬ 
clockwise  from  Ej,  and  Gj  completes  the  right-handed  system.  Ej  and  Fj  are  in 

the  Earth's  equatorial  plane,  and  Gr  is  positive  toward  the  North  Pole. 

* 
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MEAN  VERNAL 
EQUINOX 


Figure  4-1 


The  inertial  coordinate  frame  is  rotated  into  the  ECEF  system  as  follows 


E 

cos  6  sin  0  0 

EI 

F 

= 

-sin  0  cos.  9  0 

FI 

G 

R. 

0  0  1 

GI 

(31) 


I 


The  subscript  R  on  the  left-hand  side  of  equation  ( 31 ) i  indicates  that  the  EFG 
(ECEF)  coordinate  system  rotates  relative  to  the  inertial  (I)  system. 


The  reverse  transformation  is 


EI 

cos 

0 

-sin 

9 

0 

E 

FI 

sin 

9 
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GI 

I 

0 
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G 

The  transformation  of  inertial  velocities  to  ECEti 


(32) 


velocities  is 


The  reverse  transformation  is 


• 

EI 

cos 

9 

-sin  0 

0 

•  • 

E  -  8E 

• 

•  • 

FI 

V  sin 

9 

cos  9 

C 

F  +  6E 

0 

I 

o 

1 

• 

G 

The  transformation  of  inertial  accelerations  to  ECEF  accelerations  is 

E  cos  9  sin  8  0  Ej  +  29  F,  -  92Ej 

F  =  -sin  9  cos  8  ■  0  Fj  -  28  Ej  -  82Fj  (35 


The  reverse  transformation  is 


h 

COS  9 

-sin  9 

0 

••  • 

E  -  28 

• 

F  - 

82E 

fi 

= ■ 

sin  9 

cos  8 

0 

F  +  28 

• 

E  - 

92F 

i 

0 

0 

1, 

G 

The  matrix  in  equation  (31)  transforms  the  position  data  of  the  true  of  date 
inertial  frame  to  the  true  of  date  ECEF  frame  by  rotating  through  the  hour  angle 
between  the  true  vernal  equinox  and  the  Greenwich  Meridian.  This  hour  angle  is 
also  known  as  Greenwich  Apparent  Sidereal  Time  (8).  The  angle  8  can  be  computed 
from 


e  ■  cast,  ♦ <t,  -  t0hUT> 


where 


GASTq  is  the  Greenwich  Apparent  Sidereal  Time  in  radians  at  0  Universal 
Time  (t^y)  on  the  Julian  date  of  the  given  epoch, 

tj  is  the  time  of  interest,  and 
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uj  is  the  rotational  rate  of  the  Earth  measured  relative  to  the  instan- 
e 

taneously  true  vernal  equinox: 

u  =  0.7292115855  x  10‘4  .adians  per  second 
e 

cug  =  6.300388099  radians  per  day. 

These  parameters  are  illustrated  in  figure  4-2. 


Mean  Vernal  Equinox 

figure  4-2 

The  Greenwich  Apparent  Sidereal  Time  at  0 ^  UT,  expressed  as  an  angle,  repre¬ 
sents  the  angle  between  the  vernal  equinox  and  the  Greenwich  Meridian.  Assuming 
that  the  time  of  interest  (t. )  is  after  midnight,  the  Greenwich  Meridian  rotates 
through  the  angle  $  to  E'.  The  angle  9  represents  the  Greenwich  Apparent  Sidereal 
Time  at  the  time  t..  The  Greenwich  Mean  Sidereal  Time  in  \  s  at  0*1  UT  on  the 
Julian  date  epoch  day  (JDO^)  is 

GMSTq  (hours)  *  6.67170278  +  0.0657098232  (JD0£  -  2433282.5)  (38) 

MOD  (GMSTq,  24) 
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The  mean  longitude  of  the  ascending  node  of  the  Moor.'s  orbit  measured  in 
degrees  at  UT  is 

ft( Degrees)  =  372.1133  -  0.0529539  (J00E  -  2433282.5)  (39) 

Q(radians)  =  n(Degrees) 

The  equation  of  the  equinoxes  in  hours  is 

EQ  (hnurs)  =  -0.00029  sin  ft  (40) 

The  Greenwich  Apparent  Sidereal  Time  in  hours  at  0h  UT  is 

GAST  (hours)  =  GMST  +  E  (41) 

GAST  in  radians  is 
o 

GASTq  (radians)  =  GASTq  (hours)  (42) 

The  angle  9  required  for  transforming  inertial  coordinates  to  ECET  is 
G(radians)  -  GASTq  +  6.300388099  (t.  -  t0hjj)  (43) 

Mod  ( 

where 

.  ‘l  -  »,hiT  ■  (JD  -  2’4 

and  J00E  is  the  Julian  date  (JO)  at  0h 
Julian  date  of  the  time  of  interest), 
of  the  days  and  fractions  of  days  from 
Noon  (=12^UT).  The  JD  for  t.  may  be  c 


,  2tt) 


106)  -  (J00E  -  2.4  x  106) 


IT  on  the  day  of  the  epoch  (JD  is  he 
The  Julian  date  is  a  continuous  count 
1  January  4713  B.C.,  Greenwich  Mean 
imputed  u s  follows: 


(44) 


+  UT/24  -  0.5  sign  (TOOK  +  M  -  190002.5) 

+  0.5 

where 

K  is  the  year  (1801  £  K  £  2099), 

M  is  the  month  (1  £  M  <_  12), 

I  is  the  day  of  the  month  (1  <_  I  <  31), 

UT  is  the  Universal  Time  in  hours, 

<  >  is  the  integer  function,  and 
sign  is  the  sign  function. 

For  example,  1978,  January  1 ,  0h  UT  =  2443509.5  JD,  and 

1877,  August  11,  /h  30m  UT  ='2406842.8125  JD 

5.  Spherical  Inertial  Coordinate  System  (r,'  a,  5,  Vj,  Aj,  Yj,  8j) 

A  spherical  system  is  assumed  with' the  origin  the  same  as  the  EC  I  (Fj,  Fj, 
Gj)  coordinate  system.  This  coordinate  system  is  right-handed  with  Ej  positive 
toward  the  mean  vernal  equinox  of  date,  Fj  90°  counterclockwise  from  Ej,  and  Gj 
coincident  with  the  Earth's  rotational  axis  and  positive  toward  the  North  Pole 
(see  figure  5-1 ) , 
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VERNAL  EQUINOX 

Figure  5-1 

A  coordinate  system  designated  Up  Vp  Wp  is  also  shown  in  figure  5-1. 
This  coordinate  system  is  right-handed  with  Uj  pointing  east,  Vj  pointing 
north,  and  Wj  pointing  along  the  projection  of  the  geocentric  radius  r.  The 
{ U - V ) j  plane  is  normal  to  the  geocentric  radius.,  The  following  definitions 
apply  to  the  parameters  used  in  conjunction  with  the  Spherical  Inertial  Coordi¬ 
nate  System: 

r  is  the  geocentric  range  from  the  center  of  the  Earth  to  the  point 
(origin  of  the  (UVW)j  coordinate ' system) . 

5  is  the  declination,  positive  north.  The  dec 1 i nation  is  the  angle 
between  the  equatorial  plane  and'  the  range  vector. 

a  is  the  right  ascension  measured  positive  east  from  the  vernal 
equinox  of  date. 

Aj  is  the  azimuth  of  the  velocity  vector  projected  on,  the  (U-V)j 
plane. 

Yj  is  the  elevation  of  the  velocity  vector  above  the  (’J-V)j  plane. 

Bj  is  the  angle  between  the  range  vector  and  the  velocity  vector. 

If  the  position  vector  components  of  a  point  in  space  are  given  in  the  EC  I 
system,  then  r,  a,  5  may  be  found  from  thp  following: 


<EI 


cp 
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(45) 


or 


•  •  * 

Uj  =  -Ej  sin  a  +  Fj  cos  a 

Vj  *  -Ej  sin  6  cos  a  -  Fj  sin  5  sin  a  +  Gj  cos  5  (51) 

•  »  •  • 

Wj  »  Ej  cos  5  cos  a  f  Fj  cos  5  sin  a  +  Gj  sin  5 

The  transformed  velocities  may  be  used  to  find  the  total  velocity  (vj)  as 
well  as  Aj,  Yj .  3j  from  the  following  equations  (see  figure  5-2): 


Yj  *  (Uf.+  V j  +  Wf)S 


(Double  Arctan) 

(If  Aj  <  0,  then  Aj  =  Aj  +  2ir) 


(53) 


(54) 

(55) 


Figure  5-2 


If,  on  the  otner  hand,  v^,  Aj,  and  Bj  are  known,  then 

Uj  =  Vj  sin  Bj  sin 

V j  =  Vj  sin  Bj  cos  Aj  (56) 

W,  =  -v j  cos  Bj 
or 

U j  =  y j  cos  Yj  sin  Aj 

Vj  =  Vj  cos  Yj  cos  Aj  (57) 

Wj  =  Vj  sin  Yj 

In  this  case,  the  geocentric  inertial  velocities  may  be  obtained  with  the 
following  matrix  transformations: 


(58) 


Satellite  Orbital  Elements 


At  any  instant,  the  position  and  motion  of  an  artificial  satellite,  in 
earth-centered  inertially  fixed  coordinates,  can  be  described'  by  the  rectangular 
components  of  position  (Ej,  Fj,  Gj)  and  velocity  (E. ,  Fj,  Gj).  The  satellite's 

position  at  a  given  time  can  also  be  described  by  six  elements  of  the  Keplerian 
ellipse  (n,  uj,  i ,  a,  e,  M).  Three  of  these  elements  specify  the  spaclal  posi¬ 
tion  of  the  orbital  plane,  two  give  the  size  and  shape  of  the  orbit,  and  the 
sixth  relates  orbital  position  to  time.  These  elements  are  defined  as  follows: 


.1  is  the  right  ascension  of  the  ascending  node, 
ui  is  the  argument  of  perigee. 

i  is  the  incl ination  of  the  plane  of  the  orbit  with  respect  to  the  equa¬ 
torial  plane. 


a  is  the  semi-major  axis  of  the  elliptical  orbit. 
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e  is  the  eccentricity  of  the  elliptical  orbit. 

M  is  the  mean  anomaly, 
c  is  the  eccentric  anomaly. 

E|»  Fj,  Gj,  « ;*e  the  inertial  geocentric  coordinates  of  the  satellite. 

•  •  • 

•Ej,  F j ,  Gj,  are  the  inertial  geocentric  velocities  components  of  the 
satellite. 

Figure  6-1  indicates  the  coordinate  system  discussed  in  this  section. 
Figure  6-2  shows  the  relationship  of  the  position  of  a  satellite  in  its  ellip¬ 
tical  orbit  and  its  projection  onto  an  auxiliary  circle  (dotted  line).  The 
symbols  X^,  and  Yu  are  orbital  plane  coordinates,  e  is  the  eccentric  anomaly 
that  is  measured  in  the  orbital  plane  from  perigee  to  S’,  f  is  the  true  anomaly 
(not  to  be  confused  with  flattening),  and  p  is  the  radius  from  the  center  of 
the  Earth  to  the  satellite. 


Figure  6-1 


The  orbital  velocities  are 


i  =  -sin  £  \  (u 
Auj  ~  |  (1  -e  cos  e)  j 

*  m  ( (l-e£}^_cose  I  /jiVs 
oi  j  (1-e  cos  e)  |  \ay 


(61) 


Z  =  0 

a' 


where 


e  is 


y  is  the  gravitational  constant  and 

e  is  the  eccentric  anomaly,  which  must  be  obtained  by 
a  solution  of 

M  =  e  -  e  sin  e  (62) 

The  equation  is  usually  solved  by  iteration..  The  first  approximation  of 

*  M  +  e  sin  M  +  ^  e2  sin  2M 


Then 


M.j  =  e-j  -  e  sin 


AM-j  =  M  -  M1 


(63) 


AM1  • 

^E1  *  ( 1-e  cos  G|) 

Add  the  Ae^  to  G|  to  give  z^<  Repeat  until  Mn  *  M. 

The  transformation  from  orbital  coordinates  to  inertial  geocentric  coor¬ 
dinates  is 
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(64) 
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where 

P-I1  =  cos  to  cos  ft  -sin  to  sin  ft  cos  i, 

P12  =  -sin  to  cos  ft  -cos  to  sin  ft  cos  i, 

P.J2  =  sin  ft  sin  i , 

P2-|  ■  cos  to  sin  ft  +  sin  to  cos  ft  cos  i, 

P22  =  "S1n  w  sin  ^  +  cos  “  cos  ^  cos  ^ » 

P23  =  ~C0S  sin  1 » 

P31  *  sin  ft  sin  i , 

P32  =  cos  w  sin  i ,  and 

P33  -  cos  i. 

Calling  the  above  matrix  the  P  matrix,  the  transformation  from  orbital  veloCi 
ties  to  inertial  geocentric  velocities  is 


Assuming  that  the  position  and  velocity  vectors  exist  in  the  ECI  system, 
the  orbital  elements  (ft,  oj,  i,  a,  e,  M)  can  be  determined.  The  semi-major  axis 
of  the  elliptical  orbit  is 


a  ,  — HT - 

2u-r  Vj2 


(67) 
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where 


r  *  (E|  +  Fj  +  Gj)  2  and. 


VI  *  ^1  +  Ff  +  ) 2 


(68) 

(69) 


Next,  calculate  the  angular  momentum  vector  L: 


along  with 


L  -  r  x  v. 


*x  =  FIGI  •  GIFI 


\  =  GIEI  •  EIGI 


*2  =  eifi  -  FIEI 


(70) 


(71) 


The  parameter  z  defines  a  unit  vector  in  the  direction  of  L  with  components 

(72) 


where 


lx/l,  £/ L  and  yi. 


L  8  (**  +  +  *£>** 


Rotating  the  position  coordinates  through  ft  and  i  produces 
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Solving  for  G‘  gives 

G'  •  (sin  i  sin  ft)  Ej  ♦  (-sin  i  cos  ft)  Fj  f  (cos  i)  Gj  (75) 
i  qu.it  ing  the  like'  expressions  for  ?'  and  G'  gives 

<  f\-  -  sin  i  sin  ft 


*y/L  =  -sin  i  cos  ft 
»  ./ L  =  cos  i 


Solving  for  ;:  from  the  above  expressions  gives 


(76) 


(Double  Arctan) 

(If  ii  <  0,  then  Q  ■  ft  ♦  2n) 


(77) 


(78) 


The  eccentricity  (e)  can  be  calculated  from  standard  formulas  Involving  the 
semi -major  axi,  and  flattening  or  the  semi -major  and  semi-minor  axes  of  the 
ellipse.  The  eccentric  anomaly  (r)  then  follows  is 


e  cos  •  * 


MI. 

e  sin  c  -  (EjEj  ♦  FjFj  >  Gj^J/tua)**  . 

-  (m")  •  «-•*  [hs-j] 


c  •  tan' 


(79) 

(80) 

(81) 


The  term  r  1$  the  magnitude  of  the  radial  component  of  velocity.  Often 
this  quantity  Is  not  readily  available;  however,  the  radius  vector  (?)  and  the 
total  Inertial  velocity  vector  (vjj  are  generally  available.  When  this  Is  the 
case,  the  term  r  r  can  be  replaced  by  ?*Vj  (dot  product.)  In  equation  (81).  In 
any  event,  the  term  r  should  not  be  co.if-js  d  with  vf  .  Finally,  compute  the  mean 
anomaly  as 

H  •  c  *  e  sin  c.  (82) 
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if 


To  determine  the  argument  of  perigee  (u)  find  the  angle  between  the  line 
of  nodes  and  the  satellite  position  measured  in  the  orbital  plane  9. 

From  equation  (74) 

E'  =  £j  cos  Q  +  Fj  sin  Q  (83.) 

F*  =  -Ej  cos  i  sin  fi  +  Fj  cos  i  cos  fi  +  Gj  sin  i 
n  /p‘\  (Double  Arctan) 

o  -  tan  |If  9<0>  then  e  =  0  +  2tt)  (84) 


Now  find  the  true  anomaly  (f),  the  angle  between  perigee  and  the  satellite 
position: 


f  =  2  tan' 


1  (M 


‘"2  sin  T 


(Double  Arctan) 


f  =  tan' 


-l  /  L(ElEi  +  FiEi  +  GI^l)  \  (Double  Arctan) 


L2  -  yr 


(if  f  <  0,  then  f  -  f  +  2tt)  (86) 


Now  finally 


co  =  0  -  f  (If  u)  .<  0,  then  w  s  w  +  Ztr) 


(87) 


It  is  advisable  to  exercise  due  care  in  determining  the  quadrants  of  the 
angular  quantities  uj,  0,  and  f.  A  way  of  making  this  determination  is  to  com¬ 
pute  both  the,  sine  and  cosine  value  for  a  given  angle. 


Geodetic  Coordinate  SystemJA, 


The  geodetic  coordinate  system  i?  related  to  the  geocentric  coordinate 
system  as  shown  in  figure  7-1.  The  geodetic  latitude  (4>)  is  the  angle  between 
the  geodetic  Equator  and  the  geodetic  line  that  is  normal  to  the  surface  of  the 
ellipsoid  at  any  point.  The  geodetic  latitude  is  positive  in  the  northern  half 
of  the  hemisphere.  The  geodetic  longitude  {X)  of  a  point  is  measured  positive 
to  the  east  of  the  Greenwich  Meridian.  The  angle  is  measured  in  the  E-F  plane 
between  the  zero  meridian  plane  and  the  geodetic  meridian  plane  of  the  point. 
Both  planes  contain  the  minor  axis  of  the  reference  ellipsoid. 


Greenwich 

Meridian 


Figure  7-1  , 

The  geocentric  coordinates  of  a  point  in  space  may  be  obtained  from  the 
geodetic  coordinates  as  follows: 

E  =  (N  +  H)  cos  <P  cos  X 

F  *  (il  +.  H)  cos  <£  sin  X  (88) 

G  B  [N  (1-e2)  +  H]  sin  <t>,  1 

where  N.  is  the  radius  of  curvature  in  the  prime  vertical 


In  this  relation,  a  and  e  are  the  semi-major  axis  and  eccentricity,  respec¬ 
tively,  for  the  reference  ellipsoid. 


wwn» 
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IT  the  geocentric  coordinates  are  known,  the  reverse  information  is 


A  =  tan 


-l  /£]  {Double  Arctan) 

l  E j  (If  A  <  0,  then  A  =  A  +  2r) 


(90) 


Latitude  is  obtained  by  Newton's  iterative  method  as  follows: 


Vl  =  * J  '  f  J  ( t 


-l 


(91) 


where 


ae^ 


.  (E2  *  f2)^  , 


fJ  =  <1  -e-  sin2  J,)la  '  cos  *j“  '  Si" 


and  (92) 


e4a  sin  ^  cos  4>j  (E2  +  F2)'5  sin  4>j  G  cos 


j  (1-e2  sin2  4ij) 


3/2 


cos2  ip. 


sir *  4». 

J 


(93) 


For  j  =  1 


<b .  =  tan 
3 


-i 


(1-e2)  D 


(94) 


where  D  =  (E2  +  F2)'2  . 


-l 


The  solution  is  complete  when  -fj  ^  is  <  10" 12  radians.  Then 

J 


H  *  f  |(D-D$)2  +  (G  -  G$)2 


(95) 


where 


D$  «  N  cos  $  and 


(96) 


G$  »  N  (1-e2)  sin 


To  resolve  the  ambiguity  in  sign,  compute 


B  =  0s  (D-0S)  +  G$  (G-Gs) 


(97) 
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(3)  ECEF 

Spheri cal 


(2)  ECEF  Geocentric 
(Master) 


s  * 


(5)  EC  I 

Spherical 

r,  a,  4 


(4)  EC  I 

Inertial 

EI>  FI»  GI 

Fj,  Gj 


Figure  8-1 


(6)  Satellite 
Oribtal  j 

n»  '(*>,  i ,  a,  e,  M 
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VIII.  RCC  RECOMMENDATIONS  AND  PROCEDURES 

The  Range  Commanders  Council  adopts  and  recommends  fcr  interrange  use  the 
Geodetic  Datum  and  Global  Coordinate  System  described  belcw. 

a.  For  the  exchange  of  data  among  ranges  on  missions  where  interrange  action 
is  required,  employ  the  Department  of  Defense  World  Geodetic  System  of  1972,  with 
the  following  ellipsoid  definitions: 

WGS  72  ELLIPSOID 


STANDARD  ERROR 


PARAMETERS 

NOTATION  ' 

MAGNITUDE 

(68.27?) 

1* 

Gravitational 

Cohstant  , 

GM 

398600.5  km3/s2. 

■  ±0.4 

Second  Degree* 

Zonal 

^2,0 

-484.1605  x  10'6 

->  * 

Angular  Velocity*- 

ue 

0.7292115147  x  ID'1* 
rad/s 

*0.1  x  10"1 

o 

Angular  Velocity' 

u,e 

0.7292115855  x  ID"1* 
rad/s 

Sen  'ajor  Axis* 

a 

6378135m 

±5  m 

Flattening* 

f  • 

1/298.26 

±0.6  x  10' 7 

,  4 

Gravitational 

Constant  ■ 

GM'  . 

398600.8  km3/s2 

±0.4 

Major  eccentricity 

e 

0.08181331066 

b.  Transformation  constants  .or  shifting  from  three  regional  datums  to  the 
WGS-72  Datum  are  given  below.  The  unclassified  regional  datum  Shifts  listed  are 
taken  from  reference  1 . 


| 

* 

TfZ 

/.v 


r.V 


‘‘Defining  'parameters  of  *33  72  Ellipsoid  (see  reference  2 ) . 

Excluding  atmosphere  (for  use  with  geodetic  computations  involving  the  .normal 
potential). 

2IAU  system  of  Astronomical  Constants  (relative  to  fixed  stars). 

3 Relative  to  the  instantaneous  true  equinox. 

“ Contains  the  mass  of  the  Earth's  atmosphere  (for  use  with  satellite  and 
space  studies). 
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rwevtou*  »aq« 

"7* 

If  HANK 

JS 

tfy  .S  FORMAT  ION 

North  American 
Datum  (NAD)  27 
Area  to  WGS-72 

European  Datum  (ED) 
( International ) 

Area  to  WGS-72 

Tokyo  Datum  (TD) 
(Bessel ) 

Area  to  WGS-72 


i 


l...  .  . 


*1  ( m) 

AF(m) 

AG(in) 

.'-a  (m) 

Li  x  10 

-22 

157 

176 

-71.4 

-0.37295850 

-84 

-103 

-127 

-253 

-0.14223913 

-140 

516 

673 

737.845. 

0.10006272 
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3*  oarflUiath  7^’  An*  and  HN  using  aN»  fN*  FN*  and  S  as  shown  in  secti°n  VII 

pai  ayf  apn  /  • 

The  origins  and  ellipsoids  for  a  few  datums  are 


NAD  1927 

ED 

TD 

,  Meade's  Ranch.  KS 

Potsdem 

Tokyo 

♦  39*13'26.686MN 

52*22*51. 45"N 

35*39*17. 51HN 

x  261*27*29. 494"E 

13*03*  58. 93*’E 

139*44*  40. 90"E 

H  599.4M+0.3 

Ellipsoid  Clarke  1866 

International  1924 

Bessel  1841 

The  ellipsoid  parameters  of  Interest  are 

Ellipsoid 

a  (meters) 

f 

Clarke  1866 

6378206.4 

1/294.9786982 

Fischer  (Mercury)  1960 

6378166 

1/298.3 

Kaula  1961 

6378165 

1/298.3 

Bessel  1841 

6377397.155 

1/299.1528128 

MGS  72 

6378135 

l/298i26 

International  (Hayford)  1924 

6378388 

1/297.0 

MGS  84 

6378137 

1/298.257223563 
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VIII.  RCC  RECOMMENDATIONS  AND  PROCEDURES 

The  Range  Commanders  Council  adopts  and  recommends  for  Interrange  use  the 
Geodetic  Datum  and  Global  Coordinate  System  described  below. 

a.  For  the  exchange  of  data  among  ranges  on  missions  where  interrange 
action  is  required,  employ  the  Department  of  Defense  World  Geodetic  System  1984, 
with  the  following  ellipsoid  definitions: 


WGS 

84  ELLIPSOID 

PARAMETERS 

NOTATION 

MAGNITUDE 

STANDARD  ERROR 
(68.27X) 

Gravi tational** 

Constant 

GM 

398600.5  km3/s2 

+0.06 

Second  Degree* 

Zonal 

^2,0 

-484.16685  x  10"6 

+1.30  x  ir9 

Angular  Veloc1ty2* 

ue 

0.7292115  x  10*4 
rad/s 

+0.1500  x  10rJl 

Angular  Veloc1ty3* 

we 

0.72921158553  x  10-4 

+4.3(10-15  T  )  rad/s 
u 

+0.1500  x  10*11 

Gravitational  * 

Constant 

GM' 

398600.15  kmW 

+0.6 

Seml-uejor  Axis* 

a 

6378137  m 

+2  n 

Flattening  (Elllptlclty) 

f 

1/298.257223563 

(0.00335281066474) 

First  eccentricity 

e 

0.0818191908426 

e2 

0.00669437999013 

*0efln1ng  parameters  of  WGS  84  Ellipsoid  {see  reference  1). 

^Contains  the  mass  of  the  Earth's  atmosphere  {for  use  with  satellite  and  space 
studies). 

?IAG  adopted  value  for  the  true  angular  velocity  of  the  earth. 

•^Relative  to  the  instantaneous  true  equinox;  T  »  Julian  Centuries  from  Epoch 
J2000.  * 

^Excluding  atmosphere  (for  use  with  geodetic  computations  Involving  the  normal 
potential). 
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b.  Transforma tion  constants  for  shifting  from  throe  regional  datums  to  the 
WGS  84  Datum  are  given  below.  The  unclassified  regional  datum  shifts  listed  are 
taken  from  reference  1. 


TRANSFORMATION  E(m)  F(m)  G(m)  ,  a(m)  f  x  10 


North  American*  -8  160  176  -69.4  -0.37264639 

Datum  (NAD)  27 
Area  to  WGS-84 

European  Datum  (ED)  -87  -98  -121  -251  -0.14192702 

1950  ( International ) 

Area  to  WGS-84 

Tokyo  Datum  (TD)  ,  -128  481  664  739.845.  0.10037483 

(Bessel) 

Area  to  WGS-84 


*Heah  Value  (CONUS) 


Revised  June  1989 


IX.  BIBLIOGRAPHY 

1.  World  Geodetic  System  1984,  DMA  TR  8390.2,  Defense  Mapping  Agency, 
Washington,  D.C.,  30  September  1987. 

2.  Supplement  to  Department  of  Defense  World  Geodetic  System  1984  DMA 
Technical  Report:  Part  I  -  Methods,  Techniques,  and  Data  Used  in  wfis  84 
Development;  DMA  TR  5350. 2-A,  Defense  Mapping  Agency,  Washington,  D.C., 

1  December  1987., 

3.  Supplement  to  Department  of  Defense  World  Geodetic  System  1984  Technical 
Report:  Part  II  -  Parameters,  Formulas  and  Graphics  for  t.ie  Practical 
Application  otwGS  84;  DMA  Tft  8350.2-8,  Defense  Mapping  Agency, 
Washington,  D.C.,  1  December  1987. 

4.  Coordinate  Systems  and  Map  Projections,  White  Sands  Missile  Range  Data 
Systems  Manual,  September  1979,  Data  Sciences  Division,  National  Range 
Operations  Directorate,  White  Sands  Missile  Range,  NM 

5.  Bomford,  G.  GEODESY,  Oxford  University  Press,  London,  1971. 

6.  Ewing,  C.  E.,  and  M.  M.  Mitchell.  Introduction  to  Geodesy,  American 
Elevlser,  New  York,  1970. 

7.  Greenhood,  D.  Mapping,  The  University  of  Chicago  Press,  Chicago,  1964. 

8.  Haiskanen,  W.  A.,  and  H.  Moritz.  Physical  Geodesy,  W.  H.  Freeman  and 
Co.,  San  Francisco,  1967. 

9.  Hedman,  E.  L.,  A  High  Accuracy  Relationship  Between  Geocentric  Cartesian 
Coordinates  ar>d~Geodetic  Latitude  and  Altitude.  Journal  of  Spacecraft  ” 
and  Rockets,  Volume  7,  July-August  ld/0^ 

10.  Kaula,  W.  M.  Theory  of  Satellite  Geodesy,  Blalsdell  Publishing  Co., 
London,  1966. 

11.  Mallng,  0.  H.  Coordinate  Systems  and  Map  Projections.  George  Phillip  and 
Son,  London,  19/3. 

12.  Mellulsh,  R.  K.  An  Introduction  to  the ^Mathematics  of  Map  Projections. 
Cambridge  University  Press,  London,  193j. 

13.  Mueller,  I.  I.  Spherical  and  Practical  Astronomy,  Frederick  Ungar 
Publishing  Co.,  New  York,,  1969. 

14.  Olson,  0.  K.,  Optimal  Coordinate  Transformation  Formulas,  Pacific 
Missile  Range,  Point  Hcgu,'  dA,  Tft-3450-1-77.  January  1977. 

15.  Rlchardus,  P.,  and  R.  X.  Adler.  Map  Projections,  North  Holland 

Publishing  Company,  Amsterdam,  19727  _ 

lo.  Smart,  W,  A.  Spherical  Astronomy,  Cambridge  University  Press,  London, 
1965. 


